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Abstract. Inverse scattering problem is discussed for the Maxwell's equations. A reduction 
of the Maxwell's system to a new Fredholm second-kind integral equation with a scalar weakly 
singular kernel is given for electromagnetic (EM) wave scattering. This equation allows one 
to derive a formula for the scattering amplitude in which only a scalar function is present. If 
this function is small (an assumption that validates a Born- type approximation), then formulas 
for the solution to the inverse problem are obtained from the scattering data: the complex 
permittivity e'(x) in a bounded region D C K 3 is found from the scattering amplitude A(f3, a, k) 
known for a fixed k = oj^/eoMo > and all /3, a £ S 2 , where S 2 is the unit sphere in R 3 , eo and 
fio are constant permittivity and magnetic permeability in the exterior region D' = R 3 \D. The 
novel points in this paper include: 

i) A reduction of the inverse problem for vector EM waves to a vector integral equation with 
scalar kernel without any symmetry assumptions on the scatterer, 

ii) A derivation of the scalar integral equation of the first kind for solving the inverse scattering 
problem, 

and 

iii) Presenting formulas for solving this scalar integral equation. The problem of solving this 
integral equation is an ill-posed one. A method for a stable solution of this problem is given. 

Keywords and phrases: Electromagnetic waves, Maxwell's equations, wave scattering, in- 
verse scattering 
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1. Introduction 

There is a large literature on electromagnetic wave (EM) scattering ([1]- [3]), to name a few books. 
In [4] (see also [5]-[10] and references therein), wave scattering theory is developed for small bodies of 
arbitrary shapes. The idea of this paper is to reduce the Maxwell's system to one Fredholm-type integral 
equation for the vector E of cletric field and use this equation for solving inverse scattering problem. 
Let e'(x) = e(x) + icr ( x > ; where a(x) is the conductivity and e(x) is the permittivity. It is assumed that 
e(x) = 6q in D 1 = R 3 \ D, i.e., outside of an arbitrary large but finite domain D, a{x) = outside D, and 
the magnetic permeability [x = fig in R 3 , where eo and no are constants. We assume that e'(x) is twice 
continuously differentiable function. 
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The inverse scattering problem consists in finding complex permittivity e'(x) = e(x) + lcr ^ from the 
scattering amplitude A(j3, a, k), known at a fixed wavenumber k > and all directions a of the incident 
plane wave and all directions j3 of the scattered waves. In Section 2 this integral equation (2.15) is derived. 
In Section 3 equation (2.15) is used for a study of inverse scattering problem. The original ideas in this 
paper include equation (2.15) and inversion procedure in Section 3, see equations (3.8)-(3.13). 

The novel points in this paper include: 

i) A reduction of the inverse problem for vector EM waves to a vector integral equation with a scalar 
kernel without any symmetry assumptions on the scatterer, 

ii) A derivation of a scalar integral equation for solving the inverse scattering problem, and in presenting 
formulas for solving this integral equation. The problem of solving this integral equation is an ill-posed 
one. A method for a stable solution of this problem is given. 

The results are formulated in Lemmas 2.1 and 2.2 in Section 2 and Theorems 1,2 in Section 3. 

2. EM wave scattering by a body 

Let D be a finite body filled with a material that has complex permittivity e'(x) and constant magnetic 
permeability fig. The Maxwell's equations are: 

VxE = wfj,H, V x H = -iud(x)E inK 3 , (2.1) 

where u > is the frequency, /i = /io is the constant in R 3 magnetic permeability, e'(x) = eo > in D' , 
eo is a constant, 

e'{x) = e{x)+i-^-, a(x)>0, 

e'(x) ^ OVx £ R 3 , e'(x) G C 2 (R 3 ) is a twice continuously differcntiablc function, e(x) = eo m D' , and the 
conductivity a{x) = in D' . From (2.1) one gets 

VxVx£=f 2 (i)£, H = .* , (2.2) 

iujfi 

K 2 (x) := Lu 2 e'(x)fi. (2.3) 

We seek the solution of the equation 

V x V x E = K 2 (x)E (2.4) 

satisfying the radiation condition: 

E(x)=E (x)+v, (2.5) 

where Eq(x) is the plane wave, 

Eq(x) =£e lka - x , k = -, (2.6) 
c 

c = is the wave velocity in the homogeneous medium outside D, e = eo = const is the dielectric 

parameter in the outside region D', fj, = fj,Q = const in R 3 , a G S 2 is the incident direction of the plane 
wave, £ ■ a = 0, £ is a constant vector, and the scattered field v satisfies the radiation condition 

dv 1 

ikv — o(-), r = \x\ — s- oo, (2-7) 



dr 

uniformly in directions /3 := —. 

If E is found, then the pair {E,H}, where 



r 
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solves our scattering problem. 

Our goal is to derive a Fredholm second-kind integral equation for E. We follow closely in this 
derivation the presentation in [8]. Let us rewrite equation (2.4) as 

-AE + V(V • E) - k 2 E -p{x)E = 0, (2.8) 

where 

p{x) := K 2 (x) - k 2 , p{x) = in D' . (2.9) 

Note that 

JmK 2 (x) = knp(x) > 0. 

This inequality will be used in the proof of Claim 3 below (see formula (2.21)). 
It follows from (2.4) that 

= V • (K 2 (x)E) = VK 2 {x) ■ E + K 2 (x)V ■ E. (2.10) 
From (2.9) and (2.10) one gets 

—AE - k 2 E - p(x)E - V{q(x) ■ E) = 0, (2.11) 

where 

VK 2 (x) Ve'(x) . 

q{x) : = low = W q{x)= 1 ( } 

From (2.11) and (2.5) one gets 

E = Eo + J g(x, y) (p(y)E(y) + V y (q(y) ■ E(y))j dy 7 (2.13) 

where 

e ik\x-y\ 

g(x,y):=— (2.14) 

Att\x — y\ 

If q vanishes at the boundary of D, i.e., VA' 2 (i) = Vp(x) vanishes at the boundary of D, then, after an 
integration by parts in the last term of (2.13), and taking into account that 

^x9(x,y) = -W y g(x,y), 

one gets an equivalent equation: 

E = E + [ g(x, y)p{y)E{y)dy + V x f g(x, y)q(y) ■ E(y)dy := E + TE. (2.15) 

J D J D 

Since we have assumed K 2 (x) ^- and K 2 (x) £ C 2 (M 3 ), it follows that equation (2.15) is a Fredholm 
equation of the second kind, because the integral operator T in (2.15) is compact in H 1 (D), where H 1 (D) 
is the usual Sobolcv space. 

Indeed, the operator 

Bf:= f g(x,y)f(y)dy 

JD 

acts from L 2 (D) into H 2 (D), and the operator 

B l E:=p{y)E + V(q{y)-E) 

acts from H 1 {D) into L 2 {D). Thus, T acts from H 1 (D) into H 2 (D) and is, therefore, compact by the 
embedding theorem. 

Let us summarize in Lemma 2.1 the result we have proved: 
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Lemma 2.1. If K 2 (x) G C 2 (R 3 ), K 2 (x) ^ 0, K 2 (x) = k 2 > in the exterior domain D' := K 3 \ D, 
then the operatorT in (2.15) is compact in H l (D), so that equation (2.13) is of Fredholm type in H l (D). 

Let us now formulate the second auxiliary result in Lemma 2.2. After its formulation the ideas of the 
proof are briefly described, The proof consists of proving three claims, that are formulated in the course 
of the proof. 

Lemma 2.2. Equation (2.15) is uniquely solvable in iJ 1 (I?). 

Proof. Let us prove Lemma 2.2. It is sufficient to prove that the homogeneous version of equation (2.15) 
has only the trivial solution. If E solves the homogeneous equation (2.15), then E solves equation (2.4) 
and satisfies the radiation condition (2.7). The only solution to (2.4) satisfying (2.7) is the trivial solution 
E = 0. 

Let us give details and prove the above claims. 

Claim 1. A solution to (2.15) satisfies equation (2.4), and, consequently, (2.10), that is, V -(K 2 (x)E) = 

0. 

Thus, Claim 1 states that (2.15) is equivalent to the original equation (2.4), which is not at all obvious. 

Claim 2. The only solution to (2.4), (2.7) is E = 0. 

Proof of Claim 1. If E solves (2.15) and is of the form (2.5), then 

(-A - k 2 )E = p(x)E + V(q{x) ■ E). 

Rewrite this equation as 

VxVxfi- V(V • E) - V(q(x) ■ E) = K 2 {x)E, 

or 

VxVx £ -V ( K2( ^- E K :J ) K2{X) - E ) = K 2 (x)E (2.16) 

Denote 

■ (K 2 {x)E) :=#r). 

Taking divergence of equation (2.16) one gets 

-Aip - K 2 (x)ij: = in M 3 . (2.17) 

The function i/j satisfies the radiation condition (2.7). 

Claim 3. The only solution to (2.17), which satisfies condition (2.7), is ij) = 0. 

We prove Claim 3 below. Assuming that this claim is proved, we infer that ip = 0, so 

V • (K 2 (x)E) = 0, 

and equation (2.4) holds. 
Claim 1 is proved. □ 

Let us prove Claim 3. Equation (2.17) can be written as: 

— - k 2 ip - p(x)if) = in M 3 , (2.18) 

where p(x) is defined in equation (2.9), and p(x) = in D' . It is known that if k 2 > 0, p(x) <G L 2 (D), 
lmp(x) > 0, D is a bounded domain, and ip satisfies the radiation condition (2.7) and equation (2.18), 
then ip = (see, e.g., [7]). 

For convenience of the reader, we sketch the proof. From (2.18) and its complex conjugate one derives 
the relation 

tpAip - tpAip - 2ilmp(x)\'4>\ 2 = 0. (2.19) 
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Integrate (2.19) over a ball Br centered at the origin of radius R, and use the Green's formula to get 

/ (V>^p -tp^-)ds-2i f lmp(x)\iP\ 2 dx = 0. (2.20) 
J\x\=R u r or J D 

By the radiation condition (2.7) for ip one can rewrite (2.20) as 

-2ik lim / \^\ 2 ds - 2i [ lmp{x)\^\ 2 dx = 0. (2.21) 

R -^°° J\x\=R JD 

Since Imp(i) > by our assumption, and k > 0, it follows that 

lim / \i)\ 2 dx = 0. (2.22) 

R ^°° J\x\=R 

Relation (2.22) and the equation 

(A + k 2 )ijj = \x\> Rq, (2.23) 

where Br D D, imply = for \x\ > Rq- (Sec [5, p. 25]). 

By the unique continuation principle for the solutions of the homogeneous Schrodingcr equation (2.18), 
it folows that -0 = 0. Claim 3 is proved. □ 

Let us prove Claim 2. This will complete the proof of Lemma 2.2. 

If E solves (2.4) and satisfies (2.7), the the pair {E,H}, where H is defined by the second formula in 
(2.2), solves the homogeneous Maxwell's system (2.1), satisfies the radiation condition, and {E,H} is in 
Hf oc (M 3 ) . It is known (see, e.g., [3]) that this implies E = H = 0. Lemma 2.2 is proved. □ 



From formula (2.15) assuming that the origin is inside D, one gets: 

In 



E = E + e-^y(p(y)E(y) + V(q(y) ■ E))dy[l + 0(1-)} (2.24) 



as \x\ -> oo, |fy := /?. 

An alternative representation of the scattered field is given in Section 3 in formula (3.2). 

3. Inverse scattering for the Maxwell's equations 

We assume that the inhomogeneity is described by the permittivity 

/ / * , . .o-(x) 
ui 

e(x) — eo in L>' := M 3 \ D, <r(x) — in D', and D is an arbitrary large finite domain, D C Br, where 
Br is the ball of radius R centered at the origin. The origin is inside D. The inverse scattering problem 
consists of finding e'(x) from the knowledge of the scattered field measured at large distances from D. 
This scattered field is defined by the scattering amplitude. Let us use equation (2.15) for finding an 
unknown e'(x) given the scattering amplitude for all a.,/3 and a fixed k > 0. Let us define the scattering 
amplitude A(j3, a, k) by the formula 

E-E = - — A(/3,a,k)+o(-), r := \x\ -> oo, /3:=~. (3.1) 
r r r 

From equations (2.15) and (3.1) one gets: 

A<J), a, k) = i- J e-^-yp{y)E{y)dy + ^ J e^^qiy) ■ E(y)dy. (3.2) 
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Denote by [E, H) = E x H the cross product of two vectors E and H , and by (E, H) = E ■ H their dot 
product. Then equation (3.2) implies 

\p,A<fi,a,k)\ = i- f D e- ih ^p{y)%E{y)]dy. (3.3) 

If one assumes that \p(y)\ -C k 2 , and Eq = £e x , then one may replace E(y) in equation (3.3) by £e lka ' x , 
and obtain in this approximation, which is similar to the Born approximation in quantum mechanics, the 
following relation: 

4*\J3,A(p,a,k)] = \P,e\ I e ik ^>yp(y)dy. (3.4) 

JD 

A novel and practically attractive feature of equation (3.4) consists of a possibility to reduce the inverse 
problem of finding an unknown permittivity in a Maxwell's system of vector equations to solving a scalar 
equation (3.7) below. Let us assume without loss of generality that \£\ = 1. Then £]\ = sin 6*. where 
9 G [0,7r] is the angle between vectors j3 and £, 



,2 , 



sm 



and (3.4) implies 



4n([(3,A((3,a,k)},[/3,£})sm- 2 e:= f = / e^-^p^dy. (3.5) 

JD 

The right-hand side of (3.5) does not depend on 9. Thus, its left-hand side, /, remains bounded for 
9 G [0,tt]. The function 

f = f(/3,a,k)=f(k(a-{3)) 

is known since A(/3, a, k), /3, 9, and £ are known. 

Therefore, the inverse scattering problem of finding e'(x) from A(/3, a, k) is reduced to finding p = p(y) 
from the data f(/3,a,k). 

If p is found, then e' can be found, because 

6 {X) = —J^-^ (3 ' 6) 

where p = fio in K 3 , k 2 = u! 2 eofiQ, and k 2 , ui and po are assumed known. 
Thus, the inverse problem is reduced to solving the scalar equation 

^ a -^ v P (y)dy = ma,k), (3.7) 

where k = w^/eo/io > is fixed, a, /3 G S 2 run independently of each other through the unit sphere S 2 , 
and the function / = f(/3,a,k) is known on S 2 x S 2 . Let us summarize the result we have proved in 
Theorem 1. 

Theorem 1. If the scattering amplitude A(f3, a, k) is known for a fixed k > and all a, /3 G S 2 , then 
one calculates the left-hand side of equation (3.5) thus finding the scalar function f(fi,a.,k), and then 
calculates the function p{y) by solving equation (3.7). 

Let us now give a method for solving equation (3.7). We assume that the data / is given with some 
error: the exact / is unknown and the "noisy" data fs are given such that ||/ — fs\\ < 5, where 6 > is 
a small number. The problem is to find p = p(y) from equation (3.7), given fs- 

Such a problem has been studied in [5] . We formulate the results from [5] without proof. The reader 
is referred to the proofs to [5], pp. 259-274. The inversion formula from [5], p. 268, is: 

Pn{x) =\\ f(k(s' - s))h N (k(s' - s))e ik( > s '- s) - x dsds', (3.8) 
Js 2 Js 2 
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where TV is a positive integer, s' = a, s = /3, ds and ds' are surface elements of S 12 , 

h N (z) := |z|a 7V (z)fc 2 (32 7 r 4 )- 1 , (3.9) 

where 

a N (z) := [ 5 N (x)e- iz - x dx, (3.10) 
Jr 3 

and 

, \x\ 2 ,n ,N_,3/2,sinJ)-bcosb,2N+3 2k\x\ 
M*):=(l-4^) (4^) ( 6 3 /3 ) - 6 =^vT3' (3 ' n) 

where i? is the radius of the ball Br containing D. It is proved in [5], p. 268, that 

\\pn(x) -p(x)\\ lHBr) -> as N —> 00. (3.12) 

Suppose now that instead of the exact data / the noisy data fs are given, \\f — /d||L 2 (s 2 xs 2 ) < Then 
we compute the approximate solution pns(x) by formula (3.4) with fs(k(s' — s)) in place of f{k{s' — s)), 
and choose N = N(S) according to the recipe on p. 269 in [5]. Then N(S) — > 00 as 5 — > 0. With this 
choice of N(5) one proves, as in [5], p. 269, that 

\\pn(S)(x) - p(x)\\ l z(b r ) -> as 8^0. (3.13) 

These results were applied by the author also to inversion of incomplete tomographic data (see [11] and 
[12], pp. 259-264). 

Let us summarize the result in Theorem 2. 

Theorem 2. A stable solution to equation (3.7) is given by formulas (3.8)-(3.13). 
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